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ABSTRACT 

We present an analysis of the properties of optical counterparts of radio sources 
down to 1 mJy. Optical identifications have been obtained by matching together 
objects from the APM and FIRST surveys over the region 148.13° < RA(2000) < 
218.13°, —2.77° < Dec < 2.25°. Selecting radio sources down to 1 mJy, and adopting 
a uniform optical limit of bj = 21.5 we find 3176 have a counterpart in the APM 
catalogue, corresponding to 13 per cent of the radio sample. 

For bj < 20.5 we can divide radio sources into resolved radio galaxies and stellar- 
like objects (principally QSO). We find the population of radio galaxies to be mainly 
made by early-type galaxies with very red colours {bj — i? up to ~ 4) and a radio-to- 
optical ratio 10^ ^ f < 10*. The contribution of starbursting objects is negligible. In 
general QSO show r > 10'' and < 6j — i? < 1. On the basis of the R magnitudes, we 
estimate the sample of radio galaxies to be complete up to z ~ 0.3. We can therefore 
divide the whole sample of radio sources into a low-z and high-z population. The low- 
z one includes the objects identified as galaxies in the APM survey, the high-z one 
includes sources either identified as QSO or with no optical counterpart for bj < 20.5. 

We find that radio galaxies are strongly clustered and highly biased tracers of the 
underlying mass distribution. Models for the angular correlation function w{9) show 
good agreement with the observations if we assume a bias factor 6 ~ 2 at z ~ 0.3. 

Key words: galaxies: active - galaxies: starburst - Cosmology observations -radio 
continuum galaxies 



1 INTRODUCTION 

The advent of a new generation of large area radio sur- 
veys which sample the radio sky down to mJy levels (e.g. 
FIRST, Becker et al., 1995; NVSS, Condon et al., 1998; 
SUMSS, Bock et al., 1999) has introduced great advantages 
for both radio astronomy and cosmology. From the cosmo- 
logical point of view, the high surface density of sources al- 
lows studies of the large-scale structure of the Universe and 
its evolution up to high redshifts (2; ~ 4) and large physical 
scales (see e.g. Cress et al., 1996; Magliocchetti et al., 1998, 
1999). Furthermore, due to the low flux limits reached, these 
surveys sample populations of sources, such as star-forming 
galaxies, where the radio signal is produced by phenomena 
other than AGN activity. 

A complete investigation of the nature of radio sources 
and their properties can only be achieved with multi- 
wavelength follow-up. In particular optical identifications 
enable the acquisition of optical spectra, which can be used 
to derive their spectral type and redshift distribution (see 



e.g. Magliocchetti et al., 2000 - MA2000 hereafter). Great 
effort has recently been made to determine the photometric 
and spectroscopic properties of radio sources at mJy lev- 
els and fainter (Gruppioni et al., 1998; Georgakakis et al., 
1999; Sadler et al., 1999; MA2000; Masci et al., 2001). These 
studies have proven extremely useful in characterizing the 
populations of radio sources, but there are still significant 
uncertainties. For instance it has been established that there 
is starbursting population of radio sources, but there is still 
some debate about the overall fraction (see for instance the 
results from Sadler et al., 1999 as compared with those de- 
rived in Gruppioni et al., 1998; Georgakakis et al., 1999; 
MA2000). 

These follow-up studies however suffer from two great 
limitations: either they survey small areas of the sky and 
therefore generate relatively small samples of objects, or 
they lack completeness. The small number of sources leads 
to a low confidence level for any result drawn from statis- 
tical studies, while radio or optical incompleteness leads to 
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biases in the determination of the redshift distribution of 
radio sources at such low-flux levels. 

In this Paper we present an analysis of the properties 
and angular clustering of ~ 4100 radio sources brighter than 
1 mjy with optical counterparts brighter bj < 22. This sam- 
ple was obtained by matching together objects from the 
FIRST and APM surveys over ~ 350 square degrees near 
the celestial equator. Despite not having redshift measure- 
ments for these sources, their joint radio, photometric and 
morphological properties can be used to infer extremely use- 
ful information on their nature. Furthermore, for bj < 21.5, 
we obtain a catalogue of sources which is ~ 100% complete 
in the optical band (MA2000) and 80% complete in the ra- 
dio band (with a completeness level rising to 100% for radio 
fluxes S'i.4GHz > 3 mJy, Becker et al., 1995). 

The photometric and morphological properties of the 
optical identiflcations are then used to divide the whole sam- 
ple of radio sources into two well defined populations. The 
first population consists of low-redshift radio galaxies and 
the second one of high-z QSO and objects with no counter- 
part on the UKST plates. The homogeneity and complete- 
ness of the low-z sample allows us to study the clustering 
properties of these sources. 

Lastly, the area chosen for our analysis coincides with 
some of the fields observed in the 2df Galaxy Redshift Sur- 
vey (Maddox 1998, Colless, 1999). This sample of optical 
identiflcations of ~ 4100 radio sources therefore provides an 
excellent starting point for further wide-area spectroscopic 
follow-up. 

The layout of the paper is as follows: Section 2 briefly 
describes the two surveys and the data coming from them. 
Section 3 presents the procedure we adopted to match ra- 
dio and optical sources together, while Section 4 is devoted 
to the analysis of the photometric properties of the optical 
identiflcations. Section 5 examines the clustering properties 
of the sample and Section 6 summarizes our conclusions. 
Throughout the paper we will assume fio = 0.4, ho = 0.65, 
A = 0.6. 



2 THE CATALOGUES 
2.1 The FIRST survey 

The FIRST (Faint Images of the Radio Sky at Twenty cen- 
timetres) survey (Becker, White and Helfand, 1995) began 
in the spring of 1993 and will eventually cover some 10,000 
square degrees of the sky in the north Galactic cap and 
equatorial zones. The beam-size is 5.4 arcsec at 1.4 GHz, 
with an rms sensitivity of typically 0.15 mjy/bcam. A map 
is produced for each field and sources are detected using 
an elliptical Gaussian fitting procedure (White et al., 1997); 
the 5(T source detection limit is ~ 1 mJy. The astrometric 
reference frame of the maps is accurate to 0.05 arcsec, and 
individual sources have 90 per cent confidence error circles 
of radius < 0.5 arcsec at the 3 mJy level, and 1 arcsec at 
the survey threshold. The surface density of objects in the 
catalogue is ~ 90 per square degree, though this is reduced 
to ^ 80 per square degree if we combine multi-component 
sources (Magliocchetti et al., 1998). The depth, uniformity 
and angular extent of the survey are excellent attributes for 
investigating, amongst others, the clustering properties of 
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Figure 1. Number of galaxies (dotted line) and stellaj: objects 
(solid line) in the APM survey as a function of bj magnitude. 



faint sources. The catalogue derived from this survey hcis 
been estimated to be 95 per cent complete at 2 mJy and 80 
per cent complete at 1 mJy (Becker et al., 1995). 

We used the 5 July 2000 version of the catalogue which 
contains approximately 722,354 sources from the north and 
south Galactic caps. This is derived from the 1993 through 
2000 observations that cover nearly 7988 square degrees of 
sky, and includes most of the area 7''20'" < RA(2000) < 
17''20"\ 22.2° < Dec < 57.5° and 2l''20™ < RA(2000) < 
ghjo-" _2.8° < Dec < 2.2°. 



2.2 The APM survey 

The APM galaxy survey has been extensively described in 
Maddox et al., (1990a, 1990b, 1996). Briefiy, it is based on 
APM scans of 5.8° x 5.8° UKST plates, covering about 7000 
square degrees in the regions Dec < —20° in the southern 
Galactic cap and -17.5° < Dec < 2.5°, -7.5° < Dec < 
2.5° respectively in the equatorial sgp and ngp. A scan of a 
typical plate records about 300,000 images, with a limiting 
magnitude for image detection bj < 22. For each image, the 
measurements show a positional accuracy of ~ 1 arcsec and 
an isophotal magnitude accuracy of ~ 0.1 — 0.2 mag. 

The photometry from each plate is corrected so that it 
is consistent with the neighbouring plates, in the way de- 
scribed in Maddox et al. (1990b). The limit for uniform im- 
age detection over the whole survey is set by the shallowest 
plate, and is found to be bj = 21.5. 

Image profiles were then used to separate galaxies from 
stars; the resulting sample of ~ 3- 10^ galaxies (APM Galaxy 
survey) is ~ 90 — 95 per cent complete, and contamination 
from stellar objects is about 5-10 per cent over the mag- 
nitude range 16 < 6j < 20.5. For bj < 20.5 we then end 
up with two well defined subsamplcs made by two different 
classes of sources, galaxies and stellar objects. Fainter than 
this, the measured galaxy profiles do not greatly differ from 
stellar proflles and the very few images appear signiflcantly 
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Table 1. Number of sources on each UKST plate considered in this work; the first column is for all the sources before any magnitude 
cut applied to the APM catalogue (i.e. a non-uniform limit 21.5 < bj < 22), while the second one illustrates the number of sources with 
bj < 21.5 - the uniform completeness limit of the survey. Columns 3 and 4 are for the objects respectively identified as non-stellar and 
stellar-objects (both for bj < 20.5). 

Field Name # objs # objs with bj < 21.5 # galaxies {bj < 20.5) # stars {bj < 20.5) 



f853 


263346 


115857 


16096 


78028 


f854 


207832 


105631 


15057 


56021 


f855 


176135 


108062 


15715 


54046 


f856 


159556 


102762 


14170 


53507 


f857 


168511 


103053 


15031 


53824 


f858 


133900 


101223 


14028 


43720 


f859 


166885 


98839 


14201 


46849 


f860 


228542 


93941 


11893 


54350 


f861 


173320 


96414 


13003 


41883 


f862 


176040 


100065 


12137 


51545 


f863 


162137 


119390 


18089 


53059 


f864 


233215 


110147 


13233 


63483 


f865 


173488 


123105 


13544 


65765 


f866 


227948 


129765 


14689 


68062 


f867 


207689 


145523 


15225 


75443 



non-stellar. This resulting decrease in images classified as 
galaxies can be seen in Figure ^ where we plotted the distri- 
bution of non-stellar (dotted line) and stellar objects (solid 
line) as a function of the bj magnitude). 



2.3 Uniformity of the APM Data 

In order to test for systematic effects across the plates we 
have analysed the surface density of images as a function 
of position within each plate. We divided each plate into a 
grid of 15 X 15 Cij cells in RA and Dec and counted the 
number of objects in each cell, summing corresponding cells 
on different plates together. In the absence of systematic 
effects, the number of sources in each of the cells dj = 
S Plates '''J should be the same within Poisson errors, since 
superposing sources from different regions of the sky will 
average over any real structure in the distribution. 

The resulting cell counts are presented in Figure ^ 
where each line shows the cell count plotted against Dec in- 
dex j for fixed RA index i, and the different lines represent 
sets of cells with different RAs within each plate. The left- 
hand panel shows the number of objects within each set of 
cells Cij for all the sources detected i.e. with a non-uniform 
detection limit, varying between 21.5 ^ bj <^ 22. 

There are clearly systematic differences in the average 
number of images in different parts of a plate. In particular 
all the lines are lower at the ends that at the centres, showing 
that there are fewer objects at the edges of the each plate 
than in the centre. This is because of the geometrical vi- 
gnetting and differential desensitization of the photographic 
emulsion which increase the image detection limit near the 
plate edges. The effective transmission is quite flat close to 
the telescope axis but falls to < 70 per cent of the central 
value at the edges of the plates. The dashed and dot-dashed 
lines correspond to the sets of cells respectively placed on 
the outermost regions on the left-hand and right-hand side 
of each plate and are ~ 30% lower than the central val- 
ues. A correction has been applied to make the measured 
magnitudes uniform over the field, but this was applied af- 



ter the image detection process, and so we do not recover 
any images lost below the original detection limit at the 
edges of each plate. However, if we adopt a magnitude cut 
of bj = 21.5, this is brighter than the original detection limit 
even in the plate corners, and so we obtain a uniformly se- 
lected sample of images. The middle panel in Figure ^ shows 
the number of sources with bj < 21.5, and it is clear that the 
a catalogue of sources is uniform across the sky. The sub- 
sample of objects identified as galaxies also appears uniform, 
as shown in the right-hand panel of Figure |^. 

2.4 Optical colours 

The Maddox et al (1990) version of the APM data includes 
magnitude measurements only in the bj band, but Irwin et 
al. have subsequently scanned UKST R plates for the same 
area of sky. This more recent data can be found in the AP M- 
CAT database at jittp: //www, ast . cam, ac .uk/~ apmcat and 
includes both bj and R magnitudes for objects down to 
bj ~ 22, 7? ~ 21. Unfortunately for us, the APMCAT data 
was processed with stellar photometry in mind, and this 
leads to magnitude estimates that are significantly different 
to the galax;y magnitudes of Maddox et al. Re-processing all 
of the R data to make it consistent with the Maddox et al bj 
data would be a major undertaking, so we used a simple ap- 
proximation to estimate the galaxy R band magnitudes. We 
make the assumption that bj R = 6japmcat ~ ^apmcat 
which assumes that the difference between the stellar and 
galaxy magnitudes is the same in bj as in R for all plates. 
This is not exactly true but is probably not in error by more 
than 0.1 magnitudes. Thus we estimate the galaxy R mag- 
nitudes as i? = bj - [&,/apmcat - -RapmcatI- 



3 MATCHING PROCEDURE 

The FIRST and APM surveys only overlap in a relatively 
small region of the sky set on the equatorial plane between 
9*" 48'" < RA(2000) < 14'' 32'" and -2.77° < Dec < 2.25°. 
This corresponds to the UKST survey fields f853 to f867. 
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Figure 2. Number of objects within each of the 15 X 15 sets of cells Cij, obtained by adding together sources from corresponding cells 
on different UKST plates (see text for details), as a function of cell number increasing with Dec. Each curve is derived by letting the 
index j vary between 1 and 15 while keeping the other index i (indicating the RA) fixed. Different curves are obtained for sets of cells 
with different RAs. Dashed and dashed-dotted lines correspond to the sets of cells respectively placed on the left-hand and right-hand 
side of each plate. The panel on the left-hand side is for sources with bj < 22, the one in the middle for bj < 21.5, while the one on the 
right-hand side is for objects identified as galaxies with bj < 20.5. 



The number of optical sources observed within this area 
down to the magnitude limit of the APM survey is ~ 2.8- 10® ; 
the number of radio objects included in the FIRST survey 
down to 1 mjy within the same area is ~ 24, 000. Table hi 
shows the number of sources in the APM catalogue appear- 
ing in each field; first and second column are respectively 
obtained for bj < 22, magnitude limit of the survey, and 
bj < 21.5, completeness limit of the survey. Columns 3 and 
4 give the numbers of galaxies and stellar objects - obtained 
as in Maddox et al. (1990a) for bj < 20.5 - on each of the 
plates considered in this work. 

We identified the optical counterparts of FIRST radio 
sources by matching together sources from the radio sam- 
ple with those from the optical catalogue and simply taking 
pairs with separation less than a chosen radius as the correct 
identifications. The maximum acceptable radius depends on 
the accuracy of the measured positions. As mentioned in 
Section 2, the positional accuracy in the FIRST survey is 
estimated to be <0.5 arcsecs at the 3 mJy level, reaching 
the value of 1 arcsec only at the survey threshold (Becker, 
Helfand & White, 1995). This makes a mean positional er- 
ror across the survey of about ~ 0.5 arcsec; 0.5 arcsec is also 
the position accuracy for objects in the UKST plates down 
to bj ~ 21.5 (Maddox et ah, 1990a), so that the expected 
value for the radius to be used in the matching procedure is 
~ \/2 X 0.5^ ~ 0.7 arcsec. However, radio and optical refer- 
ence frames may also be offset with respect to each other. 
For example MA2000 find, for a subsample of sources drawn 
from the FIRST survey, a mean positional offset of about 0.8 
arcsec between optical (APM) and radio reference frames. 
Lastly it should also be noted that the centre of radio emis- 
sion is often displaced from the centre of optical emission, 
especially when the radio source is extended. These two fur- 
ther effects generate an extra error term which has to be 
added in quadrature to the positional accuracies in both ra- 



dio and optical catalogues, giving rise to a matching radius 
of the order of ~ 0.8 — 0.9 arcsecs. 

Tackling the problem from a more pragmatic point of 
view, we considered the distribution of the residuals Ax — 
2;aADio - a;opTicAL, Ay = i/radio - J/optical between the 
positions of all radio and optical pairs with separations | Ax| 
and \Ay\ less than 5 arcsecs. 

The results are presented in Figure where the bottom 
panel shows the Ax - Ay distribution with a point for each 
radio-optical pair, while the middle and top panels respec- 
tively show histograms of number of matches as a function 
of Ax and Ay offsets. The true optical identifications show 
up as a highly peaked concentration of points near zero off- 
set, with a rms value of the distributions (both in Ax and 
Ay) of about 0.7 arcsec. Note that this value is very close 
to the one formerly obtained by following a simple qualita- 
tive analysis, and encloses ~ 70 per cent of the candidate 
matches found within 5 arcsecs. It follows that a 2 arcsec 
match radius is equivalent to about 2.5a, and should then 
include ~ 97% of the true identifications. We therefore as- 
sume all the objects with positional residuals greater than 2 
arcsecs to be random coincidences, and choose the value of 
2 arcsecs as the correct matching radius. 

Columns 2 and 3 of Table ^ respectively give the total 
number of radio sources and number of sources with opti- 
cal identifications (with bj < 22) found in each field. The 
total number of sources in the area 148.13° < RA(2000) < 
218.13°, -2.77° < Dec < 2.25° with an identified optical 
counterpart within 2 arcsecs is 4075, 16.7 per cent of the 
whole radio sample. This percentage drops to 13 per cent 
(corresponding to 3176 identifications, column 4 of Table ^ 
if we include the further constraint bj < 21.5, to obtain a 
uniform sample from the APM survey. 

In order to quantify how many matches are likely to be 
random coincidences, we ran the matching code once again, 
after shifting all the radio positions 1 arcmin north. The 
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Table 2. Number of radio sources on each UKST plate considered in this work; the first column shows the total number of radio sources 
from the FIRST survey, while the second and third ones illustrate the number of sources with an optical identifications on the plates, 
respectively for 6j < 22 - magnitude limit of the survey - and bj < 21.5 - completeness limit of the survey. Columns 4 and 5 are for 
radio objects respectively identified as galaxies and stellar images (both for bj < 20.5) 

Field Name # radio # identifications # identifications {bj < 21.5) # galaxies {bj < 20.5) # stellar (bj < 20.5) 



£353 


1584 


304 


215 


105 


47 
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1 Kno 
IDUz 
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f855 


1531 


234 


183 


105 


29 


f856 


1578 


252 


208 


89 


53 


1857 


1589 


255 


214 


106 


45 


f858 


1579 


234 


204 


90 


41 


f859 


1690 


295 


234 


123 


39 


f860 


1588 


278 


196 


80 


55 


f861 


1576 


275 


199 


94 


36 


1862 


1582 


268 


191 


83 


49 


1863 


1695 


258 


225 


130 


42 


1864 


1730 


287 


193 


78 


45 


1865 


1613 


272 


223 


98 


64 


1866 


1788 


307 


249 


115 


44 
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1674 
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228 


98 


45 
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Figure 3. Lower panel: distribution of the residuals Ax = 
3:radio - a;oPTlCAL, Aj/ = j/radio - ^optical between radio 
and optical positions. Middle and top panels: number of matches 
as a function of Ax and Ay offsets. 



Figure 4. Lower panel: distribution of the residuals Aa: = 
a^^RADlO - a;oPTlCAL, At/' = I/radio - ^optical between ra- 
dio and optical positions, where a vertical shift of 1 arcmin has 
been added to all the radio positions. Middle and top panels: 
number ol matches as a function of Ax and Ay' offsets. 



results are shown in Figure The positional residuals are 
uniformly distributed on the Aa; Ay' plane as expected for 
random coincidences. Furthermore, we find only 225 sources 
with offsets within our chosen matching radius of 2 arcsecs, 
corresponding to ~ 5 per cent of the true matches. Con- 
tamination from spurious coincidences at this level is not 
a problem in our analysis. Note that if we had chosen a 
matching radius as large as 5 arcsecs, the percentage of ran- 
dom matches would have been 1200/5866 ~ 20 per cent, 
which would have required a more sophisticated treatment of 
contamination. These numbers are in good agreement with 
the results obtained by evaluating the expected number of 
chance coincidences: for 68.7 radio sources per square de- 



gree, one in fact respectively obtains 182 and 1139 random 
matches by considering radii of 2 and 5 arcsecs over the 350 
square degrees of analysed area. 



4 OPTICAL PROPERTIES OF RADIO 
SOURCES 

We divided the radio sources with optical identifications into 
galaxies and stellar objects based on the optical classifica- 
tions described in Section 2. Figure ^ shows the number 
of galaxies (dotted line) and stellar sources (solid line) in 
the radio subsample with optical identifications (hereafter 
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Figure 5. Number of galaxies (dotted line) and stellar-like ob- 
jects (solid line) in the radio subsample with optical identifica- 
tions as a function of bj magnitude. 
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Figure 6. Radio flux S vs. bj magnitude. The dashed lines are 
the loci of constant radio-to-optical ratio ig (see text for details). 
Red points indicate stellar-like sources (i.e. QSO), while the black 
ones are for objects identifled as galaxies. 



called the id-sample) as a function of bj magnitude. As ex- 
pected from the analysis performed in Section 2, the number 
of sources identified as galaxies rises with increasing magni- 
tude up to bj ~ 20.5, where the optical images become too 
faint to allow the detection of any extended structure. Cor- 
respondingly, the number of stellar objects increases beyond 
its real value for bj > 20.5. Again bj — 20.5 is seen as the 
faintest magnitude for which we can discriminate between 
extended images (radio galaxies) and stellar images (mainly 
QSO). The stellar identifications also include a very small 
number of galactic stars - less than 0.1 per cent of the ob- 
jects in the FIRST survey are expected to be associated with 
stars, see Helfand et al., 1997). Columns 5 and 6 of Table ^ 
indicate the number of non-stellar and stellar sources found 
on each of the UKST plates in our analysis. 

A joint analysis of the radio and photometric proper- 
ties of the id-sample can provide us with extremely useful 
insights on the nature of the radio sources under exami- 
nation. The bj and R magnitudes are plotted against radio 
fluxes in Figures |^, where black dots indicate galaxies and 
red dots indicate stellar objects. Note that beyond bj ~ 20.5 
(corresponding to i? ~ 19.5) the distinction between galax- 
ies and stars becomes very difficult, therefore all the sources 
appear in the catalogue as stellar. The following analysis will 
therefore only apply to magnitudes bj < 20.5. 

The radio-to-optical ratio is defined £LS — S X 
]^Q{m-i2.5)/2.5^ where S is the radio flux in mjy and m is 
the apparent magnitude (either R or bj depending on the 
specific case. We recall here that a source is considered 
radio-loud if tb ^ 10 (see e.g. Urry & Padovani 1995). It 
has been shown (see for instance Magliocchetti, Celotti & 
Danese, 2001) that for star-bursting galaxies, where the star- 
formation activity is so intense that it hides the signatures 
of any central BH, tb ^ 100. Thus the small number of 
sources below the r_B=100 line in Figure ^ shows that there 
are < 100 star-forming galaxies in the FIRST sample down 
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Figure 7. Radio flux S vs. R magnitude. The dashed lines are 
the loci of constant radio-to-optical ratio r/j (see text for details). 
Red points indicate stellar-like sources (i.e. QSO), while the black 
ones are for objects identifled as galaxies. 



to 1 mJy. This is further confirmation of the result that star- 
bursting sources show up only at radio fiuxes below the mJy 
level (see for instance, Gruppioni et al. 1998, Georgakakis et 
al. 1999 and MA2000). 

The region of the S-bj plane between r ~ 100 and r 
~ lO** is mainly dominated by sources identified as (radio) 
galaxies, and this population constitutes the majority of the 
id-sample. For r > lO** and bj < 20.5 stellar identifications, 
most likely to be high-z QSO, dominate the sample. Fainter 
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Figure 8. Colour bj — R vs flux for all the sources in the id-sample 
with bj < 20.5. Red points indicate stellar-like sources (i.e. QSO), 
while the black ones are for objects identified as galaxies. 



plots the bj — R colour of sources in the id-sample for ob- 
jects with bj < 20.5. The majority of the id-sample consists 
of galaxies with very red colours, up to foj — 7? ~ 4. This 
conclusion is actually made even stronger as in our analysis 
we might have missed very red sources which had bj magni- 
tudes too faint to be included in the original catalogue. The 
id-sample was drawn from a catalogue selected on bj mag- 
nitude; the red magnitudes were added later as explained in 
Section 2. Note the clear distinction in Figure M between ra- 
dio galaxies, which appear mostly with bj — R ^ 1 and QSO, 
which dominate the region < 6j — i? < 1. This division 
becomes striking for radio fluxes S > 10, where the two dis- 
tributions take the shape of a tuning fork. Again we can see 
that star-forming galaxies showing blue (i.e. for b,j ~ R'^ 0) 
colours and low (S < 10) radio-fluxes constitute only a few 
tens of objects in the lower part of the plot. 

Also when we consider the number of sources per unit 
of radio flux (see Figure ^ we see that stellar and non-stellar 
ids show different behaviours. The number of galaxies (the 
dashed line), shows a steeper slope than the number of QSO 
(the solid line) . This difference reflects the fact that galaxies 
tend to have lower radio fluxes, while QSO dominate the 
region 5'i.4GHz > 20 mJy. 
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Figure 9. Number of sources per flux interval in the id-sample. 
The dotted line represents the whole sample, while the dashed 
and solid lines respectively are for galaxies and stellar-like sources 
(with bj < 20.5) 



than bj — 20.5 this effect is mainly an artefact of the data 
because the optical images are too faint to resolve galaxies. 

Figure ^ illustrates the flux-magnitude diagram in the 
R band. Apart from the features already seen in Figure ^ 
we note an interesting shift of the population of radio galax- 
ies towards brighter red magnitudes. This shift is especially 
significant for S > 10 mJy and is not observed for the stellar 
identifications, highlighting the fact that radio-galaxies are 
typically red, early-type objects. 

This effect can be seen more clearly in Figure ^, which 



5 SPATIAL DISTRIBUTION AND 

CLUSTERING PROPERTIES OF THE 
SAMPLE 

The acquisition of optical identifications for radio sources is 
not only useful in providing optical information, it also gives 
the opportunity to make a distinction between low-redshift 
sources and more distant ones. 

Passive radio galaxies are known to be reliable stan- 
dard candles (see e.g. Hine & Longair, 1979; Rixon, Wall & 
Benn, 1991). The mean absolute magnitude in the red band 
is Mr ~ -23 (see e.g. MA2000), and there is little scat- 
ter about the mean (MA2000 estimate AM_r ~ 0.25 at the 
1(7 level). Thus measurements of the apparent magnitude R 
for these sources can provide us with an estimate of their 
redshifts according to the relation 



R — Mr = — 5 + 51ogiQdL(pc) 



(1) 



where di the luminosity distance. The id-sample is complete 
to J? ~ 19.5, and from this formula we expect that objects 
brighter than R ~ 19.5 should be found in the redshift range 
< 2 < 0.3. 

As we saw in the previous Section, the optically extended 
id-sample is mostly made up of early-type galaxies, with 
star-bursting objects contributing only a small fraction. So 
we can apply equation (^ to the whole subsample of radio 
galaxies, and deduce that it contains radio-sources that are 
mostly closer than z ~ 0.3. The spread in the R — z relation 
means that the upper redshift limit will be a gradual cutoff; 
see also MA2000. 

The other sources that appear in the id-sample are 
mainly QSO with a very small contamination from nearby 
stars. Unfortunately, we have no simple way to distinguish 
between these two classes of objects, or to estimate the QSO 
redshift. However, apart from rare exceptions, we expect 
QSO to be all placed at high redshifts, at least beyond z 
> 0.3. 
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Figure 10. Projected distribution of radio sources with optical identifications: bottom panel for the whole sample {bj < 21.5), middle 
panel for the sources identified as galaxies (bj < 20.5) and top panel for those identified as stellar (fej < 20.5). 



The optical distinction between point-like and extended 
sources as discussed in the previous Section, allows us to di- 
vide the id-sample in a low-z subsample - comprising mostly 
of early-type galaxies and roughly complete to 2 ~ 0.3 - and 
a sample of QSO mostly found at redshifts well beyond 0.3. 

Figure hd shows the angular distribution of these 
sources projected on the sky. The bottom panel represents 
the distribution of all the objects in the id-sample with 
bj < 21.5, limit for completeness and uniformity of the op- 
tical catalogue, while the middle one includes radio sources 
identified as galaxies and the top one is for point-like sources 
(both in the case of bj < 20.5). 

We can quantify the level of clustering in each sample 
by means of the two-point correlation function. Ideally one 
would like to obtain the spatial correlation function ^(r) but, 
since we do not have measured redshifts for any sources in 
the id-sample, we have to deal with its angular counterpart 
w{d). We recall here that the angular two-point correlation 
function measures the excess probability, with respect to a 
random Poisson distribution, of finding two sources in the 
solid angles 5Q2 separated by an angle 9. It is defined 
as 

5P ^11^5^15^2 [1 + w (9)] (2) 

where n is the mean number density of objects in the cata- 
logue under consideration. 

We calculated w{6) using the estimator (Hamilton, 
1993) 

^^M^ -1' (3) 

where the number of data-data, random-random and data- 
random pairs separated by an angle 9 are denoted by DD, 



RR and DR respectively. 

Random catalogues were generated with a spatial distri- 
bution modulated by both the APM and FIRST coverage 
maps, so that the instrumental window functions do not 
affect the measured clustering. We estimate the errors on 
the 1^(9) measurements by assuming the the pair counts fol- 
low Poisson statistics. This underestimates the uncertain- 
ties, but we do not expect this bias to be a large factor. 
However, the errors in w{6) at each point are correlated, 
so it is easy to overestimate the significance of any features 
apparent in w. 

In Figure |l^ we show the results for w{9) of the whole 
id-sample with bj < 21.5. As already stated, error-bars are 
given by Poisson estimates for the catalogue under consider- 
ation. The clustering signal is significantly greater than zero 
at all angular scales; if we assume a power-law form for w{9), 
w{9) = A 9^-^ , we find, via a fit to the data, A ~ 0.018 
and 7 ~ 1.3. 

However, as we have already seen, the id-sample is a hybrid 
mixture of low- and high-z sources, which means that this 
measurement is not particularly relevant. More interesting 
information can be derived from the analysis of the cluster- 
ing signal produced by the two different classes of objects, 
namely stellar and extended sources. 

We then estimated w{9) for the stellar-like sources (535 
objects with b,j < 20.5), and the resulting 1 + 'w{9) is shown 
in Figure |l^. The clustering signal is dominated by large 
errors due to the small number of objects. Nevertheless the 
signal is positive at almost all angular scales up to 6 ~ 5° . 

The most interesting result is represented by the angu- 
lar correlation function for the 1494 radio galaxies in the id- 
sample. As illustrated by figure in this case the signal is 
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Figure 11. The angular correlation function 1 + w{8) for the 
whole sample of radio sources with an optical identification within 
bj = 21.5. Error-bars are obtained as Poissonian estimates for 
the catalogue under consideration. The dashed-dotted line shows 
the best fit to the data and is obtained for a functional form 
w{e) = Ae'^-'i, where A ~ 0.018 and 7 ~ 1.3. 
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Figure 12. The angular correlation function l+w{8) for the sub- 
sample of radio sources (with bj < 20.5) showing an optical com- 
pact structure. Error-bars are as in Figure 



quite strong, showing once again that radio sources are more 
clustered than normal galaxies (see also Peacock & Nichol- 
son, 1991; Cress et al., 1996; Loan, Wall & Lahav, 1997; 
Magliocchetti et al., 1998: Magliocchetti et al., 1999). In 
fact, if we parameterize w{9) as a power-law: w{9) — AO^^'' , 
we find - via a fit to the data - A ^ 0.03 and 7 ~ 2.1, 
with a value for the slope close to those obtained for sam- 
ples of early-type galaxies only, known to be more strongly 
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Figure 13. The angular correlation function l+w{8) for the sub- 
sample of radio galaxies (with bj < 20.5). Error-bars are as in 
Figure The dashed line represents predictions from clustering 
evolution models (see text for details). 



clustered than other galaxy populations (see e.g. Maddox et 
al., 1990c and Loveday et al., 1995). Note that this result 
also agrees with the observational evidence for AGN-fuelled 
radio sources to be hosted by elliptical/early- type galaxies 
(see e.g. McLure et al., 1999). 

Since this last measurement has been obtained for a ho- 
mogeneous sample of sources with a fairly well determined 
redshift distribution, we can predict the clustering for a 
range of models to be directly compared with the data. The 
first step in our model is to calculate the mass-mass corre- 
lation function 5m(r, z). We use the analytical form derived 
using the Peacock & Dodds (1996) approach, extended to z 
7^ by Matarrese et al. (1997). This starts with a primor- 
dial CDM power spectrum (COBE normalized in our case; 
Bunn & White, 1997) and writes the non-linear part of the 
power-spectrum Pjvi(fcjvL) as a function of its linear part 
PhikL), where fci and kNL are the linear and non-linear 
wave-number respectively (see Magliocchetti et al. (200Qa) 
for a complete derivation of the various quantities). 

The predicted angular correlation function Wth{9) can 
then be obtained from ^m{r, z) using the relativistic Limber 
equation (Peebles, 1980) 



j;--N{zr^b{zrdzj;^^^{r,z) du 
[r-Niz)dzf 



(4) 



where x is the comoving coordinate, u is related to the 
spatial distance r between two sources separated by an 
angle 9 (in the small angle approximation) via r ~ 

(y? /F'^ + x^d^Y^^ {F is the curvature correction), N{z) is 
the redshift distribution for the sources under consideration, 
Zmax is the maximum redshift of the sample and b{z) is the 
bias function which determines the way radio sources trace 
the mass distribution. Note that, by imposing a Zmax we 
have assumed that the sample completeness is a step func- 



© 0000 RAS, MNRAS 000, 000-000 



10 M. Magliocchetti, S. J. Maddox 



tion, which immediately drops to zero beyond the redshift 
where the sample starts to become incomplete (0.3 in our 
case). We use the N{z) obtained by MA2000 (Fig. 13). This 
is based on redshift measurements for 76 radio sources se- 
lected from both the FIRST survey and the Phoenix sur- 
vey (Hopkins et al., 1998; Georgakakis et al., 1999) surveys. 
These samples were selected at a limit of 1 mjy and include 
redshift measurements up to z ~ 0.5. 

A little more attention has to be devoted to the issue 
of bias and its evolution. Of the various models presented in 
literature, the so-called merging model (Mo & White, 1996; 
Matarrese et al., 1997; Moscardini et al., 1998; Martini & 
Weinberg, 2001) has the most appealing physical motiva- 
tion; dark matter halos simply undergo dissipative collapse 
(i.e. merging). However this model has been shown to not 
fit the clustering data at low-redshifts (see e.g. Baugh et al., 
1999; Magliocchetti et al., 1999; Magliocchetti et al., 2000a 
to mention just a few). The reason for the discrepancy may 
be that the merging rate of objects is very high for z > 1 
but becomes negligible for z < 0.5. 

A better bias model to describe the data in this situation 
is the so-called galaxy conserving model (Fry, 1996). This 
model assumes that galaxies form at some epoch and then 
the clustering pattern simply evolves as the galaxies follow 
the continuity equation without losing their identity. It can 
be shown (Nusser & Davis, 1994; Fry, 1996) that bias for 
such galaxies evolves as 

^W = l + ^, (5) 

where D{z) is the linear growth rate for clustering, and bo 
is the bias at the present epoch, bo — o-s,RG/o's,m, where 
(78 is the rms fluctuation amplitude in a sphere of radius 8 
h~^Mpc. We calculate osmg assuming a correlation radius 
ro ~ 10 h~^Mpc for radio galaxies (Peacock & Nicholson, 
1991). We calculate ag.m using a COBE-normahzed ACDM 
model with = 0.4, ho = 0.65, A = 0.6. The resulting 
prediction for Wth(0) is presented in Figure |l^. The match 
with the data is extremely good, especially when one con- 
siders that there are no free parameters in the model. This 
gives us confidence on our understanding of the processes 
regulating clustering in general and in particular the one of 
radio sources, at least in the nearby universe. 



6 CONCLUSIONS 

We have presented an analysis of the properties of optical 
counterparts of radio sources at the 1 mJy level. The opti- 
cal identifications have been obtained by matching together 
objects from the FIRST and APM surveys over the region 
of the sky 148.13° < RA(2000) < 218.13°, -2.77° < Dec < 
2.25°. We found 4075 radio objects with fluxes greater than 
1 mJy to be identified in the APM catalogue for bj < 22 
within a matching radius of 2 arcsecs (id-sample). This cor- 
responds to 16.7 per cent of the original radio sample. A cut 
for bj = 21.5 has subsequently been applied to ensure a uni- 
form optical sample: 3176 identifications are found brighter 
than this magnitude limit. 

Such a wide area sample allows detailed statistical stud- 
ies on the photometry of radio sources. Furthermore, with 
the application of the methods introduced in Maddox et al.. 



(1990a), we managed to divide radio sources with optical 
identifications and bj < 20.5 into radio galaxies and stellar 
objects (mainly QSO with some contamination from nearby 
stars) . We find the population of radio galaxies to be almost 
exclusively dominated by early- type galaxies, with radio-to- 
optical ratios in the blue band between 10^ and 10"* and very 
red colours (up to bj~ R ^ 4). Starbursting objects make up 
only a negligible fraction of the sample. QSO tend to have 
r_B ^ 10^ and preferentially lie in the region < 6j — ii < 1 
of the colour-radio flux diagram at all radio fluxes. 

Since passive radio galaxies are known to be reliable 
standard candles, measurements of the apparent magnitude 
in the red band allowed us to divide the id-sample into 
two low-redshift and high-redshift catalogues, the flrst one 
- roughly complete to z ~ 0.3 - made by radio galaxies and 
the second one including high-z QSO. 

The completeness and homogeneity and defined redshift 
cut of the galaxy sample allowed a direct comparison be- 
tween models for the angular correlation function and the 
data. We find an extremely good match between observa- 
tional results and models if we treat radio galaxies as highly 
biased tracers of the underlying dark matter distribution 
and if we assume a bias evolution with look-back time of 
the form introduced by Fry (1996). This is in agreement with 
previous results obtained by Magliocchetti et al. (1999) for 
the whole sample of FIRST radio sources. 

Note that, as well as the results presented in this Paper, 
such a large sample of radio sources with optical identifica- 
tions forms an excellent basis for spectroscopic follow-up 
observations, which could, amongst others, fully determine 
the form of the redshift distribution of these objects at the 
mJy level out z < 0.3. We will tackle this issue in a future 
Paper where we will present results from the 2dF survey. 
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